Shock structures in time averaged patterns for the 
Kuramoto-Sivashinsky equation 

Hidetsugu Sakaguchi 
Department of Applied Science for Electronics and Materials, 
Interdisciplinary Graduate School of Engineering Sciences, 
Kyushu University, Kasuga, Fukuoka 816-8580, Japan 

February 8, 2008 



Abstract 

The Kuramoto-Sivashinsky equation with fixed boundary condi- 
tions is numerically studied. Shocklike structures appear in the time- 
averaged patterns for some parameter range of the boundary values. 
Effective diffusion constant is estimated from the relation of the width 
and the height of the shock structures. 
PACS numbers:05.45.-a, 47.54.+r 



An ordered structure can be obtained by time averages of disordered 
spatiotemporal patterns. Such ordered structures are experimentally found 
in Faraday waves, rotating thermal convection and electroconvection [1- 
3]. These ordered structures seem to come from the boundary effect. The 
Kuramoto-Sivashinsky equation is one of the simplest partial differential 
equations exhibiting spatiotemporal chaos. Average patterns of spatiotem- 
poral chaos in the Kuramoto-Sivashinsky equation were studied in [4]. Wc 
study average patterns of the Kuramoto-Sivashinsky equation with different 
boundary conditions. The equation in one dimension has the form 

h>t = ~h xx — h xxxx + — (h x ) , (1) 

where h = h(x,t) is a real function, x € [0, L], and the subscripts stand for 
derivatives. An equivalent equation is obtained for u = h x as 

Ut = u xx u xxxx + uu x . (2) 

Equation (2) possesses the Galilean symmetry(x — > x + vt,u — ► u + v). 
The Kuramoto-Sivashinsky equation exhibits spatiotemporal chaos and it is 
conjectured that the large-scale properties of the equation are described by 
the stochastic Burgers equation [5-8]. 

u t = Du xx + \uu x + £, (3) 

where D > 0, £(x,t) denotes the noise term, and the coefficient A of the 
nonlinear term is expected to be 1 owing to the same Galilean symmetry as 
the Kuramoto-Sivashinsky equation. The Fourier transform of £(x,t) 
satisfies 

<&WMO >=Tk 2 5 Kk ,5{t-t>). 

Characteristic solutions to the Burgers equation without the last noise term 
are the shock solutions: 

u(x) = Atanh(Kx), 

where k = XA/(2D). We have performed numerical simulations to seek for 
some shocklike structures in large-scale averaged patterns for the Kuramoto- 
Sivashinsky equation. We have used a simple real-space scheme obtained by 
replacing spatial and temporal derivatives by finite differences. The spatial 
discretization step is Sx = 0.25 and the temporal discretization step 6t = 
2.5 x 1CT 4 , and the system size L = 400. We have considered fixed boundary 
conditions, that is, 

u (0,t) = -U, u(L,t) = U, u x (0,t) = u x {L,t) = 0, 
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where ±U are the fixed values at the boundaries and U is changed as a 
parameter. The initial condition is 

u (x,0) = -U + Ux/(2l), forO<x<Z 

= U(x- L/2)/(L-2l), for/ <x<L-l 
= U + U(x- L)/(2l), for L-l<x<L 

(4) 

where I = 4 is assumed. The simulation is performed until t = ti = 1500 
and the temporal average is calculated as 

N 

u(x) = (1/JV) +nAt) 

n=l 

where ii = 500, At = 508t and N = (t 2 - «i)/Ai. 

Figure 1 displays snapshot patterns at t = 1500 and the averaged pat- 
tern. In the central region of the averaged pattern, u takes nearly and the 
mean slope of u is very small for U < 1.3. This type of averaged patterns is 
obtained for U = in [4,6]. 

A stationary shocklike pattern appears at U = 1.65. Since the solution 
is stationary in time, the averaged pattern and the snapshot pattern overlap 
in Fig. 1(c). This stationary solution can be understood from the absolute 
stability of a constant solution as follows. A constant solution u(x,t) = U 
is convectively unstable, but it becomes absolutely stable for U > 1.622. 
The stability of the constant solution u(x, t) = U is determined from the 
linear dispersion relation for the perturbation 5u(x, t) ~ Uke lkx ~ lujt around 
u(x,t) = U [9]. The linear dispersion relation is written as 

uj(k) = -Uk + i(k 2 — k A ). (5) 

The solution is absolutely stable under the condition that duj/dk = — U + 
2i(k — k 3 ) = and lmu(k) = —Ulmk + Re(k 2 — k A )<0 where k is generally 
a complex number. The transition occurs at U = U c where U c satisfies 

27U A - 68f/ c 2 -8 = 0, 

that is, U c = {(34 + V34 2 + 8 • 27)/27} 1/2 ~ 1.622. The convectively un- 
stable but absolutely stable solution is observed in our system, since the 
fixed boundary conditions are assumed and no disturbance is created at 
the boundaries. A solution which connects the two absolutely stable solu- 
tions u(x,t) = U and u(x,t) = —U appears as a shocklike structure with 
oscillating tails in Fig. 1(c). 
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Another type of average patterns appears for 1.3 < U < 1.6. Two nearly 
flat regions are observed and a shocklike structure connecting the two flat 
regions appears at U = 1.45 as shown in Fig. 1(b). The averaged value of 
u(x, t) in the flat region is not ±U and u(x, t) exhibits spatiotemporal chaos 
around the averaged value. We have assumed the form of the averaged 
pattern as Atanh(nx) and estimated the value of A and k. Figure 2 displays 
the averaged patterns and the corresponding profile u{x) = Atanh(Kx) at 
U = 1.33 and 1.45. As U is increased, A is increased and n is decreased. We 
have observed this shocklike structure in the time-averaged pattern even if 
the initial value of u(x, 0) is a random number between -0.02 and 0.02. We 
have checked up that the shocklike structure appears in the time-averaged 
pattern for a larger system with L = 3200. It is fairly robust, but we are not 
sure that the shocklike structure is maintained in an infinite size system. 

Figure 3 displays the height A of the shocklike structure as a function 
of U. The diamonds denote simulation results for L = 400 and the crosses 
denote simulation results for L = 1600. The shock amplitude A seems to be 
independent of the system size for sufficiently large L. The height A seems 
to increases from zero continuously, however, the transition is not clear, 
since the spatial fluctuation of the averaged pattern is large near U ~ 1.3. 
Figure 4 displays the effective diffusion constant D = XA/(2k) at L = 400 
(denoted by diamonds) and L = 1600 (denoted by crosses), where A = 1 is 
assumed. In our simulations A, n and D depends on U. We have plotted 
a relation of D(U) v.s. n(U) in Fig. 4. As n(U) is decreased, D(U) seems 
to increase. It might be interpreted that the effective diffusion constant 
tends to increase, as the length scale is longer. Sneppen et al. estimated 
the effective diffusion constant D and A as D ~ 10.5 and A = 2.32 from 
the analysis of the interface width < (h(x,t)— < h >) 2 > [8]. Our effective 
diffusion constant seems to approach about 10 as k — > 0. This result may be 
consistent with the result of Sneppen et al., although they consider a very 
large system and the effective coefficient A of the nonlinear term is not 1 in 
their analysis. 

In summary, we have numerically studied the Kuramoto-Sivashinsky 
equation with fixed boundary conditions. At both the boundaries, u(x, t) 
are fixed to ±U. If the U is smaller than a critical value, the averaged pat- 
tern takes nearly zero. If U is larger than another critical value, a stationary 
shock pattern with oscillating tails appears owing to the absolute stability. 
We have found a shocklike structure of the spatiotemporal chaos between 
the two critical values. We do not understand well the fairly stable structure 
in the spatiotemporal chaos, but it may be partly related to the weakness 
of the absolute instability. We are not sure whether the shocklike structures 
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are due to finite-size effects. 

In our finite-size simulations, the time averaged patterns take fairly def- 
inite nonzero amplitudes ±A. The time averaged pattern is approximated 
with the shock solution of the Burgers equation and the effective diffusion 
constant is estimated. The effective diffusion constant depends on the length 
scale of the shock width and it seems to approach about 10 for a very long 
scale. 

This work was supported in part by the Grant-in- Aid for Scientific Re- 
search No. 11837014 from the Ministry of Education, Science, Sports and 
Culture. 

References 

[1] B. J. Gluckman, C. B. Arnold, and J. P. Gollub, Phys. Rev. E 51, 1128 

(1995) . 

[2] L. Ning, Y. Hu, R. E. Ecke, and G. Ahlers, Phys. Rev. Lett. 71, 2216 

(1996) . 

[3] S. Rudroff and I. Rehberg, Phys. Rev. E 55, 2742 (1997). 

[4] V. M. Eguiluz, P. Alstr0m, E. Hernandez-Garcia, and O. Piro, Phys. 
Rev. E 59, 2822 (1999). 

[5] V. Yakhot, Phys. Rev. A 24, 642 (1981). 

[6] S. Zaleski and L. Lallemand, J. Phys. Lett. (Paris) 24, L793 (1985). 
[7] S. Zaleski, Physica D 34, 427 (1989). 

[8] K. Sneppen, J. Krug, M. H. Jensen, C. Jayaprakash, and T. Bohr, Phys. 
Rev. A 46, R7351 (1992). 

[9] R. J. Deissler, Physica D 25, 233 (1987). 



5 



Figure captions 

Fig. 1. Snapshot patterns (solid lines) u(x,t) at t = 1500 and the time 
averaged patterns (dashed lines) u{x) at U = 1.25 (a), 1.45 (b) and 1.65 (c). 
Fig. 2. Time averaged patterns (solid lines) and their fitting with u = 
Atanh(Kx) (dashed lines) at U = 1.33 (a) and 1.45 (b). 
Fig. 3. Relation of the height A v.s. U for the time averaged patterns at 
L = 400 (denoted by diamonds) and L = 1600 (denoted by crosses). 
Fig. 4. Relation of the effective diffusion constant D v.s. k for the time 
averaged patterns at L = 400 (denoted by diamonds) and L = 1600 (denoted 
by crosses). 
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